ON STRONG SOLUTIONS FOR POSITIVE DEFINITE JUMP 

DIFFUSIONS 



EBERHARD MAYERHOFER, OLIVER PFAFFEL, AND ROBERT STELZER 



Abstract. We show the existence of unique global strong solutions of a class of stochastic 
differential equations on the cone of symmetric positive definite matrices. Our result includes 
affine diffusion processes and therefore extends considerably the known statements concerning 
Wishart processes, which have recently been extensively employed in financial mathematics. 

Moreover, we consider stochastic differential equations where the diffusion coefficient is 
given by the a-th positive semidefinite power of the process itself with 0.5 < a < 1 and 
obtain existence conditions for them. In the case of a diffusion coefficient which is linear in 
the process we likewise get a positive definite analogue of the univariate GARCH diffusions. 



1. Introduction 

A result of the general theory for affine Markov processes on the cone St of symmetric 



positive semidefinite matrices developed in 13j is that for a d x d matrix-valued standard 



Brownian motion B, d x d matrices Q and (3, a symmetric constant drift b, and a positive 
linear drift T : St — > St, weak global solutions exist to the stochastic differential equation 
(SDE) 

(1.1) dX t = ^/T t dB t Q + Q T dBj^/x~ t + (Xt/3 + f3 T X t + T(X t ) + b)dt, 

X = x £ S+, 

whenever b — (d — 1)Q T Q G St. Above \fX denotes the unique positive semidefinite square 
root of a matrix X S St. For T = solutions to the SDE (|l.ip are called Wishart processes 
and their existence has been considered in detail in the fundamental paper by Marie-France 
Bru pj. Further pro babilistic investigations on properties of Wishart processes have been 
carried out in 1^, 20, 26], for instance, and references therein. 



In the present paper, we focus on the existence of global strong solutions of (jl.ip and 
generalisations of it including jumps and more general diffusion coefficients. Because of the 
non-Lipschitz diffusion at the boundary of the cone, this problem is a quite delicate one - 
a-priori it is only clear that a unique local solution of (jl.ip exists until Xt hits the boundary 
of Sj, since the SDE is locally Lipschitz in the interior of S^ . Furthermore, known results for 
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pathwise uniqueness, for instance, that of the seminal paper of Yamada and Watanabe |45l . 
Corollary 3], are essentially one-dimensional, and therefore do not apply. Hence, the present 
setting seems to be more complicated than, for instance, the canonical affine one (concerning 
diffusions on R™ x R n , |H, Lemma 8.2]). 

Positive semidefinite matrix valued processes are increasingly used in finance, particularly 
for stochastic modellin g of multivariate stochastic volatility phenomena in equity and fixed 
income models, see (^.Ticl 14-17, 24, 25, 27, H(J. See also [HI] and the references therein. 



Most papers mentioned use Bru's class of Wishart diffusions, as this results in multivariate 
analogues of the popular Heston stochastic volatility model and its extensions, Ornstein- 
Uhlenbeck type processes C [4011 ) giving a multivariate generalisation of the popular model of 
|3] or a combination of both ([3ll]). This motivated the research of 13] on positive semidefinite 
affine processes including all the aforementioned models and generalising the results of 211 ] . 
which covered all of these models in the univariate setting. Appropriate multivariate models 
are especially important for issues like portfolio optimisation, portfolio risk management and 
the pricing of options depending on several underlyings, which are heavily influenced by the 
dependence structure. 

Clearly £j~-valued processes model the covariances, not the correlations, which are, how- 
ever, preferable when interpreting the dependence structure. The results of the present paper 
are particularly relevant, when one wants to derive correlation dynamics (see e.g., 0, fiol]). 
because one needs to assume boundary non-attainment conditions for a rigorous derivation. 

The name "Wishart process" is, unfortunately, not always used in the same way in the 
literature. We follow the above cited applied papers in finance and call any solution to 
(jl.ip with r = "Wishart process" whereas in most of the previous probabilistic literature 
"Wishart process" also means (3 = and the "Wishart processes with drift" of [2(| are not 
even special cases of our "Wishart processes" . For V = j3 = and b = nQ T Q with n £ N one 
may also speak of a "squared Ornstein-Uhlenbeck process" . In the univariate case the name 
"Wishart process" is not used, instead one typically uses "Cox-Ingersoll-Ross process" in the 
financial and "squared Bessel process" in the probability literature. 

However, in this paper we do not limit ourselves to the analysis of Instead, as a 

special case of a considerably more general result, we consider a similar SDE allowing for a 
general (not necessarily linear) drift T and an additional jump part of finite variation. This 
implies that many Levy-driven SDEs on Sj" like the positive semidefinite Ornstein-Uhlenbeck 
(OU) type processes (see 0,[3§]) or the volatility process of a multivariate COGARCH process 
(see [43f] ) , where the existence of global strong solutions has previously been shown by path- 
wise arguments, are special cases of our setting. Thus our results allow to consider certain 
"jump diffusions" (in the sense of [12]) , viz. mixtures of such jump processes and Wishart 
diffusions, in applications. 

It should be noted that 0] also contains results on strong solutions for Wishart processes 
(see our upcoming Proposition 13. II and Remark I4.8p . however, they are derived under strong 
parametric restrictions, because her method requires an application of Girsanov's theorem. 
The latter is based on a martingale criterion, which in the matrix valued setting seems hard 
to verify. Also, the general result (with a non-vanishing linear drift) only holds until the first 
time when two of the eigenvalues of the process collide. Our approach generalises her method 
of proof for the case f5 = (vanishing linear drift) and avoids change of measure techniques. 

The most general result of our paper, Theorem I3.4|, also opens the way to use positive 
semidefinite extensions of the univariate GARCH diffusions of [36( or of so-called generalised 
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Cox- Ingersoll- Ross models (cf. e.g. [y, |22j), where the square root in the diffusion part of 
(jl.ip is replaced by the a-th positive semidefinite power with a £ [1/2, 1], in applications (see 
Corollary I3.5() . 

The remainder of the paper is structured as follows. In the subsequent section we summarise 
some notation and preliminaries. In Section [3] we state our main result, Theorem I3.4| and its 
corollaries applying to Wishart processes, matrix-variate generalised Cox-Ingersoll-Ross and 
GARCH diffusions. Moreover, we compare our results to the work of Bru which is recalled 
in Proposition 13. 1[ In the following section we gradually develop the proof of our result. Our 
method relies on a generalisation of the so-called McKean's argument, but avoids the use of 
Girsanov's theorem. In Section [4. II we thus provide a self-contained proof of a generalisation 
of McKean's argument and then deliver the proof of Theorem 13.41 in Section T4.2I We conclude 
the paper with some final remarks in Section 



2. Notation and general set-up 

We assume given an appropriate filtered probability space (f2, F, P, (J-'t)teR + ) satisfying the 
usual hypotheses (complete and right-continuous filtration) and rich enough to support all 
processes occurring. For short, we sometimes write just Q when actually referring to this 
filtered probability space. B is a d x d standard Brownian motion on ft and d € N always 
denotes the dimension. Furthermore, we use the following notation, definitions and setting: 

• R_|_ := [0, oo), Mfi is the set of real valued d x d matrices and Id is the identity matrix. 

• Sd C Md is the space of symmetric matrices, and S't C Sd is the cone of symmetric 
positive semidefinite matrices in Sd and S^ + its interior, i.e. the positive definite 
matrices. The partial order on Sd induced by the cone is denoted by ■<, and x >~ 0, 
if and only if x G <Sj~ + . We endow Sd with the scalar product (x,y) := Tr(xy), 
where Tr(A) denotes the trace of A £ Md- || • || denotes the associated norm, and 
d(x,dSj) = inf y( - ds + \\x — y\\ is the distance of x G S^ to the boundary dS^. 

• The usual tensor (Kronecker) product of two matrices A, B is denoted by A ® B and 
the vectorisation operator mapping Md to R rf2 by stacking the columns of a matrix 
A below each other is denoted by vec(^4) (see [29|, Chapter 4] for more details). 

• A function / : S^ + — > U with U being (a subset of) a normed space is called locally 
Lipschitz if ||/(x) — f(y)\\ < K(C)\\x — y\\\/x,y E C for all compacts C C U. / is 
said to have linear growth if ||/(x)|| 2 < K(l + ||x|| 2 ) Vx G S^ + . 

• An S^-valued cadlag adapted stochastic process X is called S'j'-increasing, if X t >z X s 
a.s. for all t > s > 0. Such a process is necessarily of finite variation on compacts by 
0, Lemma 5.21] and hence a semimartingale. We call it of pure jump type provided 
X t = X + J2o<s<t Ax °, where AX S = X s - X s _. 

For the necessary b ackg round on stochastic analysis we refer to one of the standard ref- 



erences like [3CJ, |41|, |42|. Moreover, we frequently employ stochastic integrals where the 
integrands or integrators are matrix- or even linear-operator valued. Thus, we briefly ex- 
plain how they have to be understood. Let (A t ) te ^+ in Md, (B t ) ts ^+ in Md be cadlag and 
adapted processes and (L t ) igK + in Md be a semimartingale (i.e. each element is a semi- 
martingale). Then we denote by Jjj A s _dL s B s _ the matrix Ct in Md which has r/-th element 
Cij t t = Ylt=i Yli=i Jo Aik,s-Bij, s -dLki lS . Equivalently such an integral can be understood in 
the sense of j35l | by identifying it with the integral J Q * A s ^dL s with At being for each fixed 
t the linear operator Md —> Md, X i— > AfXBt and L being a semimartingale in the Hilbert 
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space Ma- Stochastic integrals of the form f Q f K(X s J)dJ s with J being a semimartingale in 
Md (coordinatewise or equivalently as in [35j, Section 10] where the equivalence easily follows 
from [35I Section 10.9] and by noting that on a finite dimensional Hilbert space all norms 
are equivalent) and K{x) : M d — > M d a linear operator for all x can be understood again as 
in (3BI ]. Alternatively, one can equivalently identify Md with M d2 using the vec-operator and 
K(x) with a matrix in M d i d 2 and then define the stochastic integral coordinatewise as above. 

3. Statement of the main results 
3.1. Wishart diffusions with jumps. 

In order to illustrate the context of our result and, because it is of most relevance in ap- 
plications, we discuss first the special case of Wishart diffusions with jumps. For Q G Md, 
5 > d — 1, f3 G Md and an M^-valued standard Brownian motion B, a Wishart process is the 
strong solution of the equation 

(3.1) dX t = ^/x t dB t Q + Q T dBj^/x~ t + (X t p + f3 T X t + 5Q T Q)dt, 

X = x€ S++, 

on the maximal stochastic interval [0,T X ), where T x is naturally defined as 

T x = mi{t>0 : X t £dS+}. 

That such a unique local strong solution, which does not explode before or at time T x , exists, 
follows from standard SDE theory, since all the coefficients in (|3.ip are locally Lipschitz and 
of linear growth on <Sj~ + . To be more precise, this follows by appropriately localising the 
usual results as e.g. in [Il|, Chapter V] or by variations of the proofs in jH, Chapter 3]. 
A localisation procedure adapted particularly to certain convex sets like S d is presented in 



detail in [4J, Section 6.7]. 

The following is a summary of the results 0, Theorem 2, 2' and 2"] - the to the best of 
our knowledge only known results regarding strong existence of Wishart processes: 

Proposition 3.1. Let 6 > d + 1. 

(i) If Q = Id and (3 = 0, then T x = 00. 

Suppose additionally that the d eigenvalues of x are distinct. 

(ii) If Q € Sj[ + , —(3 G S~d such that f3 and Q commute, then there exists a solution 
(Xt)teR+ °f (|3.ip until the first time t x when two of the eigenvalues of X t collide. 

(iii) If (3 = (3oId and Q = jld, where (3q, 7 G R, then T x = 00 for the solution of (X t )t<=M. + 
of (jXTD- 

Consequently, for the respective choice of parameters, there exist unique global strong S^ + - 
valued solutions of the SDE f)3. 1 1) on [0, t x ) resp. on all of [0, 00). 

The upcoming general Theorem 13.41 implies the following result for a generalisation of the 
Wishart SDE allowing for additional jumps and a non-linear drift T. 

Corollary 3.2. Let b G S d , Q G M d , f3 G M d , and let 

• J be an Sd-valued cddldg adapted process which is Sf -increasing and of pure jump 
type, 

• T : S d + — > Sf be a locally Lipschitz function of linear growth and 

• K : — > L(S d ,S d ) (the linear operators on Sd mapping S^ into S d ) be a locally 
Lipschitz function of linear growth. 
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Ifby(d+ 1)Q T Q, then the SDE 

(3.2) dX t = ^QldBtQ + Q T dBj^X^+ (X t _f3 + /3 T X t _ + T(X t _) + b)dt + K(X t ^)dJ t , 

X =xeS++, 

has a unique adapted cadlag global strong solution (Xt)t^R + on . In particular we have 
T x := wf{t > : Xt~ G dSj or X t Sj + } = inf{t > : X t _ £ dSj} = oo almost surely. 

Proof. For the term on the right hand side of the upcoming condition (|3.3p we obtain 

Tr(2/3) +Tr(r(x)x" 1 ) + Tr((6- (d + l)Q T Q)x~ 1 ) > 2Tr(/3), 

noting that x — 1 , T(x) and 6— (d+ 1)Q T Q are positive semidefinite and that is a selfdual 
cone, which implies that Tr(zy) > for any z,y G S'j - . Setting c(i) = 2Tr(/3) an application 
of Theorem 13.41 concludes. □ 

By choosing T linear and J = 0, we obtain a result for (jl.ip which considerably generalises 
Proposition 13.11 

Remark 3.3. (i) in the univariate case the condition b >z (d + 1)Q T Q is known to be 
also necessary for boundary non- attainment (see Chapter XI]). 

(ii) A possible choice for J is a matrix subordinator without drift (see fl/j, i.e. an S~f- 
increasing Levy process. By choosing F / m (|3,2p appropriately our results also 
apply to SDEs involving matrix subordinators with a non-vanishing drift. 

(iii) Setting Q = 0, T = 0, K to the identity and b equal to the drift of the matrix 
subordinator, Equation (|3.2p becomes the SDE of a positive definite O U type process, 
0> \SQ] - Likewise, it is straightforward to see that the SDE of the volatility process Y 
of the multivariate COGARCH process of f^3 / is a special case of (|3,2p . 

(iv) An OU-type process on the positive semidefinite matrices is necessarily driven by a 
Levy process of finite variation having positive semidefinite jumps only (follows by 
slightly adapting the arguments in the proof of \39l. Theorem 4-9])- This entails that a 
generalisation of the above result to a more general jump behaviour requires additional 
technical restrictions. 

3.2. The general SDE and existence result. 

The main result of this paper is the following general theorem concerning non-attainment of 
the boundary of S~f and the existence of a unique global strong solution for a generalisation 
of the SDE (jl.ip . The proof of this result is gradually developed in the next sections. 

Theorem 3.4. Let 

• F,G : IR+ x S^ + — > M d , be functions such that G T <X> F given by G T &> F(t, x) = 
(G(t,x)) T (8> F(t,x) is locally Lipschitz and of linear growth, 

• H : M + x — > Sd be locally Lipschitz and of linear growth, 

• J be an S^-valued cadlag adapted process which is S~f -increasing and of pure jump 
type, 

• and K : S^ + — > L(S^ , S'j") (the linear operators on Sj mapping St into S^ ) be a 
locally Lipschitz function of linear growth. 

Suppose that there exists a function c : IR+ — > K which is locally integrable, i.e. J Q S \ c(t)\dt < oo 
for all s £ M + , such that 

(3.3) c(t) < Tr(H(t, x)x~ 1 ) - Tr(/(i, x)x~ 1 ) Tr(g(t, x)x^) - Tr(/(i, x) X - 1 g(t, x)x~ l ) 
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for all x G Sj + and t G R+ where f(t,x) := F(t,x)F(t,x) T , #(i,a;) = G(t, x) T G(t, x). 
Then the SDE 

(3.4) dX t =F(t,X t _)dB t G(t,X^) + G(t,X t „) T dBjF(t,X t _) T 

+ H(t,X t -)dt + K(X t -)dJ t , 
X =x G Sj + , 

has a unique adapted cddldg global strong solution (Xt)teR+ 071 $ d + ■ 

In particular, we have T x := mf{t > : Xt- G dSf or Xt S d + } = inf{i > : Xt- G 
dSj} = oo almost surely. 

3.3. Positive definite extensions of generalised Cox-Ingersoll-Ross processes and 
GARCH diffusions. 

In the univariate case generalised Cox-Ingersoll-Ross (GCIR) processes given by the SDE 
dxt = (b + axt)dt + qxfdB t with b > 0,q > 0,a £ U. and a G [1/2, 1] are - as discussed in 
the introduction - of relevance in financial modelling, a = 1/2 corresponds, of course, to the 
already discussed Bessel case, whereas a = 1 gives the so-called GARCH diffusions. Given 
the popularity of the Wishart based models in nowadays finance, it seems natural to consider 
also positive semidefinite extensions of the GCIR processes. An application of our general 
theorem to the case where F(X) = X a , G(X) = Q with a G [1/2, 1] yields: 

Corollary 3.5. (i) Let a G [1/2, 1], b G S d , Q G M d , /3 G M d , and let 

• J be an S d -valued cddldg adapted process which is S~t -increasing and of pure jump 
type, 

• r : Sf — > Sf be a locally Lipschitz function of linear growth and 

• K : — > L(Sf,Sj) (the linear operators on S d mapping into S j~ ) be a locally 
Lipschitz function of linear growth. 

Suppose that for all x G 

(3.5) Tr(T(x)x- 1 + bx~ l ) > Tr^ 20 " 1 ) Tr(Q T Qx~ 1 ) + Tr(x 2a ~ 2 Q T Q). 
Then the SDE 

(3.6) dX t = X?_B t Q + Q T dBjX^ + (X t -P + /3 T X t _ + T(X t -) + b)dt + K(X t _)dJ t , 

X = x G S++, 

has a unique adapted cddldg global strong solution (Xt)teisu on 5"j" + . In particular we have 
T x := inf{t > : X t - G dSj or X t = inf{t > : Xt- G 8S+} = oo almost surely, 

(ii) Any of the following sets of conditions implies ()3.5)) : 

(a) b + r(x) h Tr(x 2a ~ 1 )Q T Q + x a ~ 1 / 2 Q T Qx a ~ 1 / 2 for all x G S+ + . 

(b) b + r(x) t. Tr(x 2a - 1 )Q T Q + \ Q T Q x 2a ~ l for all x G Sj[ + with X qTq denoting the 
largest eigenvalue ofQ T Q. 

(c) a = 1 and b + T(x) >z Tr(x)Q T Q + XqtqX for all x G . 

(d) 6^0 and T(x) £. 2 Tv(x 2a - 1 )Q T Q for all x G Sj + . 

(e) b y and T(x) t 2 (Tr(x) + d(2a - l) 2 ~ 2a ) Q T Q for all x G Sj~ + (and setting 
0° := 1 for a = 1/2). 

(f) 6^0 and T{x) y 2(Tr(x) + d)Q T Q for all x G Sj + . 

(g) a > 1/2, d = 1, T(x) > for all x G R+ and b > 0. 
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Proof. One easily calculates the right hand side of (|3.3p to be equal to Tr(2/3 + T(x)x 1 + 
bx' 1 ) - Tr(x 2a ~ 1 )Tr(Q 1 ' Qx' 1 ) - Tr(x 2a_2 Q T Q) and hence (i) follows from Theorem I3U 

Turning to the proof of (ii) using the selfduality of S'f as in the proof of Corollary 13.21 gives 
(a). Next we observe that Q T Q ^ XgTgld and, hence, x a ~ 1/2 Q T Qx a ~ 1/2 ^ AqtqX 2 " -1 . 
This gives (b) and (c) is simply the special case for a = 1. 

Since for A,B£ we have that Tr(AB) < Tv(A)Tr(B) due to the Cauchy-Schwarz 
inequality and the elementary inequality \/a + b < ^/a + \fb for all a, b G M + , we have that 
Tv(x 2a ~ 2 Q T Q) < Tr(x 2a - 1 )Tr(Q T Qx- 1 ). Hence, ([23]) is implied by Tr^^)^- 1 + bx' 1 ) > 
2Tr(x 2a_1 )Tr(Q T Qx- 1 ). Using once again the selfduality gives (d). 

Since the trace is the sum of the eigenvalues, A > A 2 "" 1 for all A > 1 and a G [1/2, 1] and 
\2a-i < \ _|_ maXAe j ^ jA 2 "^ 1 — A} for all A G [0, 1) and a G [1/2, 1], we immediately obtain 
(e) from (d), because max^[ 01 ] |A 2a_1 — A} = (2a — l) 2-2 ". In turn (f) follows from (e) 
noting that max Ag [ 01 ] {A 2a_1 — A} G [0, 1]. 

Turning to (g) we have for the right hand side of (|3.3|) in the univariate case 

£(x) = 2(3 + r(x)/x + b/x - 2Q 2 /x 2 ~ 2a . 

Now one notes that the second term is non- negative and that for b > the term b/x — 
2Q 2 /x 2 ~ 2a is bounded from below on IR + , because lim x _;.o,a;>o x~ 1 /x 2a ~ 2 = oo. Hence, The- 
orem 13.21 concludes . □ 

In the different cases of (ii) a valid choice of b and V is always obtained by taking them 
equal to the right hand side of the inequalities. It should be noted that (c) shows that in 
the positive semidefinite GARCH diffusion generalisation one can always take a linear drift. 
Likewise, (e) and (f) show that a linear drift is possible for the generalized CIR. For a = 1/2 
the case (d) is again sharp in the univariate setting, but for general dimensions it is a stronger 
condition than the one given in Corollary 13.21 

The last case (g) in particular recovers the well-known univariate result for dxt = (b + 
axt)dt + qxfdBt with b > 0, q > 0, aGM and a G [1/2,1]. In our matrix-variate case 
for a > 1/2 a result similar to the univariate one, viz. that a strictly positive constant 
drift is all that is needed to ensure boundary non-attainment, seems to be out of reach. 
When one tries to use arguments similar to (e) in general, one would need something like 
Tr(6a; _1 ) > kTr(x 2a ~ 1 ) Tv(Q T Qx -1 ) + K with some constants k > and K to ensure (|3.5p . 
However, when the process comes close to the boundary of the cone, this only means that at 
least one eigenvalue gets close to zero. Hence, Tr(6x _1 ) and Tr(Q T Qx _1 ) should then go to 
infinity at a comparable rate. However, all the other eigenvalues of x may still be arbitrarily 
large and so there is no appropriate upper bound on the term Tr(x 2Q_1 ). 

4. Proofs 

In this section we gradually prove our main result. As a priori all processes involved are 
only defined up to a stopping time, we collect first some basic definitions regarding stochastic 
processes defined on stochastic intervals following mainly [33|. 

Definition 4.1. Let A £ J 7 and let T be a stopping time. 

• A random variable X on A is a mapping A — > R which is measurable with respect to 
the a-algebra in J. 

• A family (Xt)t£R + of random variables on {t < T} is called a stochastic process on 
[0, T). If X t is {t < T} n Tt-measurable for all t G K+, then X is said to be adapted. 



8 



E. MAYERHOFER, O. PFAFFEL, AND R. STELZER 



• An adapted process M on [0, T) is called a continuous local martingale on the in- 
terval [0, T) if there exists an increasing sequence of stopping times (T n ) n6 pj and a 
sequence of continuous martingales (M^) n ^ (in the usual sense on [0,oo)) such 
that lining oo T n = T a.s. and Mt = on {t < T n }. Other local properties for 
adapted processes on [0, T) are defined likewise. 

• A semimartingale on [0, T) is the sum of a cadlag local martingale on [0, T) and an 
adapted cadlag process of locally finite variation on [0, T) . 

• For a continuous local martingale on [0,T) the quadratic variation is the KU {oo}- 
valued stochastic process [M, M] defined by 

[M, M] t = sup[M (n) , M (n) ] iATn for all t G K+. 

nGN 

4.1. McKean's argument. 

In this section we finally establish Proposition 14. 31 which generalises an argument of j34l . p. 47, 
Problem 7] concerning continuous local martingales on stochastic intervals used, for instance, 
in @,@,|33]- We keep the tradition of referring to it as McKean's argument. Since it may also 
be helpful in other situations, we state our result and its proof in detail. 

Lemma 4.2. Let M be a continuous local martingale on a stochastic interval [0, T). Then 
on {T > 0} it holds almost surely that either lim t ^r M t exists in E or that lim sup^ T M t = 
— lim vaitYT = 00. 

Proof. Combine [33I . Theorem 3.5] with analogous arguments to the proof of (42I . Chapter V, 
Proposition 1.8]. □ 

Proposition 4.3 (McKean's Argument). Let Z = (Z s ) se ^ + be an adapted cadlag IR + \{0}- 
valued stochastic process on a stochastic interval [0, To) such that Zq > a.s. and tq = inf{0 < 
s < To : Z s _ = 0}. Suppose h : M + \{0} — > M is continuous and satisfies the following: 

(i) For all t G [0, tq), we have h{Z t ) = h(Zo) + M t + Ft, where 

(a) P is an adapted cadlag process on [0, To) such that '^te[o,T AT) Pt > —00 a - s - f or 
each T G M+\{0}, 

(b) M is a continuous local martingale on [0,to) with Mq = 0, 

(ii) and lim^ h(z) = —00. 
Then tq = 00 a.s. 

Above, To = inf{0 < s < To : Z s - = 0} is not to be understood as the definition of To, 
but it means that the already defined stopping time to is also the first hitting time of Z s _ at 
zero. Since Z is only defined up to time to, one cannot take the infimum over R + . 

Proof. Since h(Zt)~ = h(Zt—) = h(Zo) + + Mt- and Pt- is a.s. bounded from below 
on compacts, we have to = inf{s > : M s _ = —00} and further To > due to the right 
continuity of Z. Assume, by contradiction, that To < 00 on a set A G J- with ¥(A) > 0. 
Hence, lim^ Mt = —00 on A and this contradicts Lemma 14.21 □ 

4.2. Proof of Theorem 3.4. 

Before we provide a proof of Theorem 13.41 we recall some elementary identities from matrix 
calculus and provide some further technical lemmata. For a differentiable function / : M^ —> 
K, we denote by V/ the usual gradient written in coordinates as (^r-)ij- 

Lemma 4.4. On S~t + , we have 
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(i) Vdet(x) = det(x)(x- 1 ) T = det(x)x _1 , 
( u ) nr det ( x ) = det(x)[(x~ 1 ) H (2;~ 1 )ij - (x" 1 )^^" 1 )^]- 



Proof. The first identity in (i) can be found in [32L Section 9.10] and the second is an im- 
iate consequence of restricting to symme 

1 (dlfci' 3 ') x ~ 1 and fi nan y the symmetry: 



mediate consequence of restricting to symmetric matrices. Now (ii) follows using q§^x 



d d ( d 

det(x) = - 5 — (det(x) (x~ 1 ) ji ) = det(x) i(x~)ik(x~)ji + (x~ )ji 



dxkiXij dxki 

= det(x) ((x~ 1 )ik(x~ 1 ) ji - (x~ 1 ) jk (x~ 1 )ii) . 



□ 



For a semimartingale X we denote by X c as usual its continuous part. All semimartingales 
in the following will have a discontinuous part of finite variation, i.e. ^o<s<t ll^^ll i s finite 
for all t £ R + . Thus we define X£ = X t — X^o<s<i ^-X s and note that the quadratic variation 
of a semimartingale is the one of its local continuous martingale part plus the sum of its 
squared jumps. 

The continuous quadratic variation of X solving (|3.4p is only influenced by the Brownian 
terms and, hence, we have a general version of 0, Equation (2.4)] which is proved just as [l|, 
Lemma 2]: 

Lemma 4.5. Consider the solution Xt of f|3.4|) on [0, T x ). Then 

dlXij^XuFt = ( FF T^ Xt _ )) . k{G r G ^ Xt _ )) , i + ( FF T (* J X t _)) iI (G T G(t, X f _)) ifc 
+ (FF T (t,X t .)) jk (G T G(t,X t .)) u + (FF T (t,X t .)) jl (G T G(t,X t -)) lk . 
Here G T G(t,x) : = G(t,x) T G(t,x) and FF T (t,x) := F(t,x)F(t,x) T to ease notation. 

Moreover, we shall need the following result where a Brownian motion on a stochastic 
interval [0,T) is defined as a continuous local martingale on [0,T) with [/3,/3]j = t. 

Lemma 4.6. Let Xt be a continuous -valued adapted cadlag stochastic process on a sto- 
chastic interval [0,T) with T being a predictable stopping time and let h : — > Then 
there exists a one- dimensional Brownian motion /3 h on [0, T) such that 

(4.1) Tr (J* h{X u ^)dB^j = J y/Tr(h(X u -) T h{X u -))d/3* 

holds on [0, T). 

Proof. We define for t G [0, T), 

oh ._ h(X u -)jj ^ R 

' t^ x h ^Tr(h{X u _Yh(X u _)) UJU 

and since the numerator equals zero, whenever the denominator vanishes, we use the conven- 
tion g = l. Clearly for each i,j and for all u £ [0, T) we have 
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which ensures that f3 h is well-defined, square-integrable and a continuous local martingale on 
[0, T) by stopping at a sequence of stopping times announcing T. Furthermore, by construc- 
tion 

p*> a = y f ^.t^lL = t 



Tr(h(X u ^) T h(X u -)) 



and therefore (3 h is a Brownian motion on [0, T). 
Finally by the very definition of j3 h , we have 

d I 

Tr(h(X t -)dB t ) = Y h(X t ..) ij dB t)ji = ^TY(h{X t _yh(X t _))d$, 

which proves identity f)4. 1|) . □ 

Finally, we state a variant of Ito's formula which we later employ. It follows easily from 
the usual versions like Theorem 3.9.1] by arguments similar to [H, Theorem 5.4] and 0, 
Proposition 3.4]. 

Lemma 4.7. Let X be an -valued semimartingale on a stochastic interval [0,T) and 
f : + — > K a twice continuously differentiate function. If X t - G for all t G [0, T) and 
12o<s<t ll^-^sll < 00 / or ^ ^ [0,T) ; £/ien f(X) is a semimartingale on [0, T) and 

f(X t )=f(X )+Tr( Vf(X s ^ T dX c s )+- / «— *w]S 

\7o / i,j,k,l=i ^ Xkl 

+ Y (f(X s )-f(X s _)). 

0<s<t 

We are now prepared to provide a proof of Theorem 13. 41 Note that to shorten our formulae 
we use in the following differential notation and not integral notation as above. 



Proof of Theorem \3.4\ Since 

vec(F{t,X t -)dB t G(t,X t -)) = (G(t, X t _) T ® (F(t, X t -)) vec(dB t ), 

it is easy to see that all coefficients of (j3.4|) are locally Lipschitz and of linear growth. Hence, 
standard SDE theory implies again the existence of a unique cadlag adapted non-explosive 
local strong solution until the first time T x = mf{t > : Xt- G dS^ or Xt G" $d + } when X 
hits the boundary or jumps out of . Hence, we have to show T x = oo. 

By the choice of K and J, all jumps have to be positive semidefinite and hence the solution 
X cannot jump out of 5*j" + . This implies that T x = inf{i > : X t - G dS^}. 

In the following, all statements are meant to hold on the stochastic interval [0, T x ). Note 
that by the right continuity of Xt, a.s. T x > 0. Moreover, we set T n = inf{i G M+ : 
d(X t ,dSj) < 1/n or \\X t \\ > n}. Then (T n ) ne pj is an increasing sequence of stopping times 
with linin^oo T n = T x , hence T x is predictable. 

We define the following processes and functions according to the notation of Proposition 
~^31 



(4.2) Z t := det(Xt), h(z) : = ln(z), r t := h(Z t ). 

Then T x = inf{i > : = — oo}. 
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By Lemma 14.41 (i) and using the abbreviation / = FF , g = G G, we obtain 

TY(V(det(AV))^ t c ) =det(X t _) 2^Tr (f(t, X t _)X;}g{t, X t _)X;})dW t 

+ Tr(H(t,X t _)X^})dt , 

with some one-dimensional Brownian motion W on [0, T x ), whose existence is guaranteed by 
Lemma 14,61 Furthermore, by Lemma 14.41 (ii) Lemma 14.51 and elementary calculations we 
have that 

■det(X t _)d[X l3 ,X kl } c t 



It. 



2 dxijdxki 



^^zl UiX^MX^h - (X^) u {X^) 3k ){f(t ) X t .) lk g{t,X t .) jl 



i,j,k,l 



+ f{t,X t .) u g(t,X t _) jk + f{t,X t _) jk g{t,X t _) a + f(t,X t .) jig {t,X t _) lk ) 

= det(X t _) (Tr(f(t,X t _)X^}g(t,X t _)Xf}) - Tr(f(t,X t _)X^) Tr(g(t, X t „)X;})) dt. 
According to Ito's formula, Lemma [4. 7\ we therefore obtain by summing up the two equations, 



dZ t =2det(X t _)y/Tr(f(t,X^)Xf}g(t,X^)X^)dW t + det(X t ) - det(X t . 
+ det(AV) Tr(H(t, X t „)X^}) + Tr(/(t, X t _)X^}g{t, X t _)X~}) 

- Tr(/(t, Xt^Xf}) Tr(g(t, X t _)X^) 
Using again Ito's formula, we have 



dr t =2JTr(f(t,X t _)X^_ 1 g(t,X t _)X^_ 1 )dW t +Hdet(X t )) - ln(det(X t _)) 



+ 



Tr(H(t, X t -.)X£) - Tr(/(t, X t -)X^g(t, X^Xf}) 



dt. 



- Tr(/(t, X t -)XrJ) Tr(g(t, X^Xf}] 
Hence, we have r± = ro + Mt + Pt, where 

M t =2 ^Tr(/( S , X a -)Xj}g(s, X s „)X;})dW s , 

Tr(H(s,X s _)X;}) - Tr(f(s,X s -)X:}g(s,X s -)X-i] 



Pt 



Tr(f(s, X S -)X-}) Tr(g(s, X^)X~}) 



ds+ (ln(det(X s )) - ln(det(X s _))) . 



0<s<t 



We infer that (M t (n) ) t > defined by 

-t 



:=2f ^Tr(f(s,X^)(X^)-ig(s,X^)(Xj^)dW s 
Jo 
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is a continuous martingale. Obviously, Mt = on {t < T n } and thus M is a continuous 

local martingale on [0, T x ). Furthermore, X s — X s _ >z for all s € [0, T) and hence det(X s ) > 



det(X;_) using 28, Corollary 4.3.3]. Therefore, we have that P t > J* c(s)ds on [0,T X ). 

Finally, by Proposition 14.31 we have that T x = oo a.s. noting that c is assumed to be locally 
integrable. □ 

Remark 4.8. Bru's method for proving her proposition \3.1\ for Wishart diffusions consists 
of the following two steps: 

(i) First assume (3 = 0. By applying the original McKean's argument twice, one derives 
that /i(det(X)) is a local martingale. This is proved separately for 5 = d + 1 and 
5 > d + 1 by choosing h(z) = ln(z) in the first case and h(z) = z d+1 ~ s in the second 
one. Therefore, the existence of a unique global strong solution on + is settled. 

(ii) One may therefore suppose that Xt is an -valued solution on [0, oo) of 

dX t = y/x t dB t Q + Q T dBj^/Y t + 5Q T Qdt, X = x £ S j 



J d ■ 

where Q E GL{d) and 5 > d + 1. Now, Girsanov's Theorem is applied which allows 
to introduce a drift by changing to an equivalent probability measure. This step gen- 
eralises a one- dimensional method by Pitman and Yor, see p. 718]. The involved 



arguments and calculations, which are not presented in detail in jj/, appear rather 



complicated and work seemingly only in the special case given in Proposition \3.1\ (ii) 



m 



The technical details of |3/ concerning strong solutions are explained in more detail in fH/. 

Our proof above circumvented the problems associated to the use of Girsanov's theorem 
by extending the approach outlined in (i). 

5. Conclusion 

In this paper we have extended the previously known sufficient boundary non-attainment 
conditions for certain Wishart processes to more general SDEs on S^ + , which include affine 
diffusions with state-independent jumps of finite variation. This allowed to infer the existence 
of strong solutions of a large class of affine matrix valued processes. Moreover, we have thus 
obtained strong existence results for SDEs which can be considered as positive semidefinite 
extensions of GARCH diffusions and generalised Cox-Ingersoll-Ross processes. 

However, this results in several open questions related to our SDE (jl.ip which will hopefully 
be addressed in future work. The following questions are beyond the scope of the present pa- 
per, since they are obviously rather non-trivial and apparently need very different techniques 
than the ones employed here. For d > 1 and the Wishart diffusions it is not clear, whether 
the condition b y (d + 1)Q T Q for the drift is a necessary non-attainability condition or not. 
Only in the case j3 = 0, T = 0, Q = Id and b = 51^ with 6 £ (d — 1, d + 1) it is known from 
(20I . Theorem 1.4] that the boundary is hit. On the other hand, one knows that in the case 
d = 1 pathwise uniqueness holds, hence there exists a strong solution for all b y (even in 



the general setting of CBI processes, see 18|, Theorem 5.1]). For d > 2, the situation seems in 
general to be rather complicated and therefore existence of global strong solutions remains an 
open problem when b )£ (d+ 1)Q T Q (and the conditions for the existence of weak solutions of 
[131 ] are satisfied). Likewise, it is a very interesting problem in the case of the GCIR processes 
with a > 1/2 whether a state dependent drift away from the boundary is really necessary 
and what happens if one has only a constant drift towards the interior of Sj". 
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Finally, we remark that our method of proof could be generalised to state-spaces D other 
than St, as long as the existence of an appropriate function g : D — > IR + is guaranteed, 
such that g~ 1 (0) = dD. For instance, similar (but simpler) arguments to the ones of the 
pro of of Theorem 13.41 yield a rigorous proof of the non-attainment condition formulated in 
[111 . Section 6] for affine jump diffusions on the canonical state space IR+ x IR n . Here one takes 
g(xi,x 2 , . . -,x m ) =xt-x 2 x m . 



References 

[1] A. Ahdida and A. Alfonsi. Exact and high order discretization schemes for 
Wishart processes and their affine extensions. Preprint, 2010. available at: 
|http : //adsabs . harvard . edu/abs/2010arXivl006 . 2281A 

[2] O. Barndorff-Nielsen and V. Perez-Abreu. Matrix Subordinators and related Upsilon 
Transformations. Theory Probab. Appl., 52:1-23, 2008. 

[3] O. E. Barndorff-Nielsen and N. Shephard. Non-Gaussian Ornstein-Uhlenbeck-based mod- 
els and some of their uses in Financial Economics (with discussion). J. R. Stat. Soc. B 
Statist. Methodol, 63:167-241, 2001. 

[4] O. E. Barndorff-Nielsen and R. Stelzer. Positive-Definite Matrix Processes of Finite 
Variation. Probab. Math. Statist, 27:3-43, 2007. 

[5] K. Bichteler. Stochastic Integration with Jumps, volume 89 of Encyclopedia of Mathe- 
matics and its Applications. Cambridge University Press, Cambridge, 2002. 

[6] M. Borkovec and C. Kliippelberg. Extremal behavior of diffusion models in finance. 
Extremes, 1:47-80, 1998. 

[7] M.-F. Bru. Wishart Processes. J. Theoret. Probab., 4:725-751, 1991. 

[8] M.-F. Bru. Diffusions of Perturbed Principal Component Analysis. J. Multivariate Anal, 
29:127-136, 1989. 

[9] A. Buraschi, A. Cieslak, and F. Trojani. Correlation Risk and the Term Structure of 

Interest Rates. Working paper, Imperial College, UK, 2007. 
[10] A. Buraschi, P. Porchia, and F. Trojani. Correlation Risk and Optimal Portfolio Choice. 

J. Finance, 65:393 - 420, 2010. 
[11] P. Cheridito, D. Filipovic, and M. Yor. Equivalent and Absolutely Continuous Measure 

Changes for Jump-Diffusion Processes. Ann. Appl. Probab., 15:1713-1732, 2005. ISSN 

1050-5164. 

[12] R. Cont and P. Tankov. Financial Modelling with Jump Processes. CRC Financial 

Mathematical Series. Chapman & Hall, London, 2004. 
[13] C. Cuchiero, D. Filipovic, E. Mayerhofer, and J. Teichmann. Affine Processes on Positive 

Semidefinite Matrices. Ann. Appl. Probab., 21:397-463, 2011. 
[14] J. Da Fonseca, M. Grasseli, and C. Tebaldi. Option Pricing when Correlations are 

Stochastic: An Analytical Framework. Review of Derivatives Research, 10:151-180, 2007. 
[15] J. Da Fonseca, M. Grasseli, and F. Ielpo. Hedging (Co)Variance Risk with Variance 

Swaps. Working paper, University St.Gallen, 2008. 
[16] J. Da Fonseca, M. Grasseli, and F. Ielpo. Estimating the Wishart Affine Stochastic 

Correlation Model using the Empirical Characteristic Function. Working paper, 2008. 
[17] J. Da Fonseca, M. Grasseli, and C. Tebaldi. A Multifactor Volatility Heston model. J. 

Quant. Finance, 8:591-604, 2008. 
[18] D. A. Dawson and Z. Li. Skew Convolution Semigroups and Affine Markov processes. 

Ann. Probab., 34:1103-1142, 2006. 



14 E. MAYERHOFER, O. PFAFFEL, AND R. STELZER 

[19] C. Donati-Martin. Large Deviations for Wishart Processes. Probab. Math. Statist., 28: 
325-343, 2008. 

[20] C. Donati-Martin, Y. Doumerc, H. Matsumoto, and M. Yor. Some Properties of the 

Wishart Processes and a Matrix Extension of the Hartman- Watson Laws. Publ. Res. 

Inst. Math. Sci., 40:1385-1412, 2004. 
[21] D. Duffie, D. Filipovic, and W. Schachermayer. Affine Processes and Applications in 

Finance. Ann. Appl. Probab., 13:984-1053, 2003. 
[22] V. Fasen, C. Kliippelberg, and A. Lindner. Extremal behavior of stochastic volatility 

models. In A. N. Shiryaev, M. R. Grossinho, P. E. Oliveira, and M. L. Esqui'vel, editors, 

Stochastic finance, pages 107-155, New York, 2006. Springer. 
[23] D. Filipovic and E. Mayerhofer. Affine Diffusion Processes: Theory and Applications. In 

Advanced Financial Modelling, volume 8 of Radon Ser. Comput. Appl. Math. Walter de 

Gruyter, Berlin, 2009. 

[24] C. Gourieroux and R. Sufana. Derivative Pricing with Wishart Multivariate Stochastic 

Volatility. J. Bus. Econom. Statist, 28:438-451, 2010. 
[25] C. Gourieroux and R. Sufana. Wishart Quadratic Term Structure. Working paper, 

CREST, CEPREMAP and University of Toronto, 2007. 
[26] P. Graczyk and L. Vostrikova. The Moments of Wishart Processes via Ito calculus. 

Theory Probab. Appl, 51:609-625, 2007. 
[27] M. Grasselli and C. Tebaldi. Solvable Affine Term Structure Models. Math. Finance, 18 

(1):135-153, 2008. 

[28] R. A. Horn and C. R. Johnson. Matrix Analysis. Cambridge University Press, Cambridge, 
1985. 

[29] R. A. Horn and C. R. Johnson. Topics in Matrix Analysis. Cambridge University Press, 
Cambridge, 1991. 

[30] J. Jacod and A. N. Shiryaev. Limit theorems for stochastic processes, volume 288 of Fun- 
damental Principles of Mathematical Sciences. Springer- Verlag, Berlin, second edition, 
2003. 

[31] M. Leippold and F. Trojani. Asset Pricing with Matrix 

Affine Jump Diffusions. Working paper, 2010. available from: 

http : //papers . ssrn. com/ sol3/papers . cfm?abstract_id=1274482, 

[32] J. R. Magnus and H. Neudecker. Matrix differential calculus with applications in statistics 
and econometrics. Wiley, Chichester, 1988. 

[33] B. Maisonneuve. Une mise au point sur les martingales locales continues definies sur un 
intervalle stochastique. In Seminaire de Probabilites, XI (Univ. Strasbourg, Strasbourg, 
1975/1976), volume 581 of Lecture Notes in Math., pages 435-445, Berlin, 1977. Springer. 

[34] H. P. McKean. Stochastic Integrals. Academic Press, 1969. 

[35] M. Metivier and J. Pellaumail. Stochastic Integration. Academic Press, New York, 1980. 
[36] D. B. Nelson. ARCH models as diffusion approximations. J. Econometrics, 45:7-38, 
1990. 

[37] J. R. Norris, L. C. G. Rogers, and D. Williams. Brownian Motions of Ellipsoids. Trans. 

Amer. Math. Soc, 294:757-765, 1986. 
[38] O. Pfaffel. Wishart Processes. Project thesis, Zentrum Mathematik, Technische Univer- 



sitat Miinchen, 2008. available from: http://www-m4.ma.tum.de 



[39] C. Pigorsch and R. Stelzer. On the Definition, Stationary Distribution and Second 
Order Structure of Positive Semi-definite Ornstein-Uhlenbeck type Processes. Bernoulli, 
15:754-773, 2009. 



ON STRONG SOLUTIONS FOR POSITIVE DEFINITE JUMP DIFFUSIONS 



15 



[40] C. Pigorsch and R. Stelzer. A Multivariate Ornstein-Uhlenbeck Type Stochastic Volatil- 
ity Model. Working paper, Technische Universitat Miinchen, 2009. available from: 
|http : 77 www-m4 . ma . turn . de. 

[41] P. E. Protter. Stochastic Integration and Differential Equations, volume 21 of Stochastic 
Modelling and Applied Probability. Springer-Verlag, Berlin, 2005. Second edition. Version 
2.1, Corrected third printing. 

[42] D. Revuz and M. Yor. Continuous Martingales and Brownian motion, volume 293 of 
Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathe- 
matical Sciences]. Springer-Verlag, Berlin, 1991. 

[43] R. Stelzer. Multivariate COGARCH(l,l) Processes. Bernoulli, 16:80-115, 2010. 

[44] R. J. Stelzer. Multivariate Continuous Time Stochastic Volatility Models Driven 
by a Levy Process. PhD thesis, Centre for Mathematical Sciences, Mu- 
nich University of Technology, Garching, Germany, 2007. available from: 
http : / / mediatum2 . ub . turn . de/ node?id=624065. 

[45] S. Watanabe and T. Yamada. On the Uniqueness of Solutions of Stochastic Differential 
Equations II. J. Math. Kyoto Univ., 11:553-563, 1971. 

Vienna Institute of Finance, University of Vienna and Vienna University of Economics and 
Business Administration, Heiligenstadterstrasse 46-48, 1190 Vienna, Austria 
E-mail address: eberhard.mayerhof er@vif .ac.at 

TUM Institute for Advanced Study & Zentrum Mathematik, Technische Universitat Munchen, 
Boltzmannstrasse 3, D-85747 Garching bei Munchen, Germany 
E-mail address: pfaffel@ma.tum.de 

Institute of Mathematical Finance, Ulm University, Helmholtzstrasse 18, D-89081 Ulm, Ger- 
many 

E-mail address: robert.stelzer@uni-ulm.de 



